Quantifying single gene copy number by measuring fluorescent probe lengths on combed genomic DNA

An approach was developed for the quantification of subtle gains and losses of genomic DNA. The approach relies on a process called molecular combing. Molecular combing consists of the extension and alignment of purified molecules of genomic DNA on a glass coverslip. It has the advantage that a large number of genomes can be combed per coverslip, which allows for a statistically adequate number of measurements to be made on the combed DNA. Consequently, a high-resolution approach to mapping and quantifying genomic alterations is possible. The approach consists of applying fluorescence hybridization to the combed DNA by using probes to identify the amplified region. Measurements then are made on the linear hybridization signals to ascertain the region's exact size. The reliability of the approach first was tested for low copy number amplifications by determining the copy number of chromosome 21 in a normal and trisomy 21 cell line. It then was tested for high copy number amplifications by quantifying the copy number of an oncogene amplified in the tumor cell line GTL-16. These results demonstrate that a wide range of amplifications can be accurately and reliably quantified. The sensitivity and resolution of the approach likewise was assessed by determining the copy number of a single allele (160 kb) alteration.

rrndb: the Ribosomal RNA Operon Copy Number Database

The Ribosomal RNA Operon Copy Number Database (rrndb) is an Internet-accessible database containing annotated information on rRNA operon copy number among prokaryotes. Gene redundancy is uncommon in prokaryotic genomes, yet the rRNA genes can vary from one to as many as 15 copies. Despite the widespread use of 16S rRNA gene sequences for identification of prokaryotes, information on the number and sequence of individual rRNA genes in a genome is not readily accessible. In an attempt to understand the evolutionary implications of rRNA operon redundancy, we have created a phylogenetically arranged report on rRNA gene copy number for a diverse collection of prokaryotic microorganisms. Each entry (organism) in the rrndb contains detailed information linked directly to external websites including the Ribosomal Database Project, GenBank, PubMed and several culture collections. Data contained in the rrndb will be valuable to researchers investigating microbial ecology and evolution using 16S rRNA gene sequences. The rrndb web site is directly accessible on the WWW at http://rrndb.cme.msu.edu.

Simultaneous A8344G heteroplasmy and mitochondrial DNA copy number quantification in Myoclonus Epilepsy and Ragged-Red Fibers (MERRF) syndrome by a multiplex Molecular Beacon based real-time fluorescence PCR

The association of a particular mitochondrial DNA (mtDNA) mutation with different clinical phenotypes is a well-known feature of mitochondrial diseases. A simple genotype–phenotype correlation has not been found between mutation load and disease expression. Tissue and intercellular mosaicism as well as mtDNA copy number are thought to be responsible for the different clinical phenotypes. As disease expression of mitochondrial tRNA mutations is mostly in postmitotic tissues, studies to elucidate disease mechanisms need to be performed on patient material. Heteroplasmy quantitation and copy number estimation using small patient biopsy samples has not been reported before, mainly due to technical restrictions. In order to resolve this problem, we have developed a robust assay that utilizes Molecular Beacons to accurately quantify heteroplasmy levels and determine mtDNA copy number in small samples carrying the A8344G tRNALys mutation. It provides the methodological basis to investigate the role of heteroplasmy and mtDNA copy number in determining the clinical phenotypes.

Comprehensive copy number and gene expression profiling of the 17q23 amplicon in human breast cancer

The biological significance of DNA amplification in cancer is thought to be due to the selection of increased expression of a single or few important genes. However, systematic surveys of the copy number and expression of all genes within an amplified region of the genome have not been performed. Here we have used a combination of molecular, genomic, and microarray technologies to identify target genes for 17q23, a common region of amplification in breast cancers with poor prognosis. Construction of a 4-Mb genomic contig made it possible to define two common regions of amplification in breast cancer cell lines. Analysis of 184 primary breast tumors by fluorescence in situ hybridization on tissue microarrays validated these results with the highest amplification frequency (12.5%) observed for the distal region. Based on GeneMap'99 information, 17 known genes and 26 expressed sequence tags were localized to the contig. Analysis of genomic sequence identified 77 additional transcripts. A comprehensive analysis of expression levels of these transcripts in six breast cancer cell lines was carried out by using complementary DNA microarrays. The expression patterns varied from one cell line to another, and several overexpressed genes were identified. Of these, RPS6KB1, MUL, APPBP2, and TRAP240 as well as one uncharacterized expressed sequence tag were located in the two common amplified regions. In summary, comprehensive analysis of the 17q23 amplicon revealed a limited number of highly expressed genes that may contribute to the more aggressive clinical course observed in breast cancer patients with 17q23-amplified tumors.

Quantitative analysis of chromosomal CGH in human breast tumors associates copy number abnormalities with p53 status and patient survival

We present a general method for rigorously identifying correlations between variations in large-scale molecular profiles and outcomes and apply it to chromosomal comparative genomic hybridization data from a set of 52 breast tumors. We identify two loci where copy number abnormalities are correlated with poor survival outcome (gain at 8q24 and loss at 9q13). We also identify a relationship between abnormalities at two loci and the mutational status of p53. Gain at 8q24 and loss at 5q15-5q21 are linked with mutant p53. The 9q and 5q losses suggest the possibility of gene products involved in breast cancer progression. The analytical techniques are general and also are applicable to the analysis of array-based expression data.

High copy number I-Ab transgenes induce production of IgE through an interluekin 4-dependent mechanism.

To better understand the role of class II major histocompatibility complex molecules in both normal and autoimmune responses, we have produced a series of I-Ab transgenic mice. One of these transgenic constructs, designated NOD.PD, has the sequence of the NOD beta chain (Abeta(g7)) except at positions 56 and 57, where Pro-Asp replaces His-Ser. Several NOD.PD transgenic lines have been produced. One line of these mice carried a very high number of copies (>50) of the NOD.PD transgene. As has been described in other mice carrying high copy numbers of I-Ab transgenes, B-cell development was abnormal. The steady state numbers of mature B cells (IgM+/IgD(hi)) in the periphery were greatly reduced in transgenic mice compared to nontransgenic littermates. Surprisingly, rather than being accompanied by a generalized hypogammaglobulinemia, this B-cell deficiency was accompanied by elevated concentrations of IgG1 and IgE in the serum. Conversely, the levels of IgG2a were reduced in transgenic mice compared to nontransgenic littermates. Because this isotype pattern was characteristic of interleukin (IL)-4-induced class-switching, we then investigated the role of IL-4 in causing the observed phenotype. We crossed the high copy number transgenic mice with an IL-4-deficient strain of mice. As expected, the elevated levels of IgE in high copy number transgenic mice were eliminated when the IL-4 gene was inactivated. However, the reduction in the number of B cells was not ameliorated. These data indicate that the primary defect caused by the transgene was to reduce the number of B cells in these mice. This reduction was accompanied by a secondary increase in IL-4 production, which drove the remaining B cells toward the production of IgGl and IgE.

Single-copy transgenic mice with chosen-site integration.

We describe a general way of introducing transgenes into the mouse germ line for comparing different sequences without the complications of variation in copy number and insertion site. The method uses homologous recombination in embryonic stem (ES) cells to generate mice having a single copy of a transgene integrated into a chosen location in the genome. To test the method, a single copy murine bcl-2 cDNA driven by either a chicken beta-actin promoter or a human beta-actin promoter has been inserted immediately 5' to the X-linked hypoxanthine phosphoribosyltransferase locus by a directly selectable homologous recombination event. The level of expression of the targeted bcl-2 transgene in ES cells is identical in independently isolated homologous recombinants having the same promoter yet varies between the different promoters. In contrast, the expression of bcl-2 transgenes having the same (chicken beta-actin) promoter varies drastically when they are independently integrated at random insertion sites. Both promoters direct broad expression of the single-copy transgene in mice derived from the respective targeted ES cells. In vitro and in vivo, the human beta-actin promoter consistently directed a higher level of transgene expression than the chicken beta-actin promoter.

Chromosome-specific molecular organization of maize (Zea mays L.) centromeric regions

A set of oat–maize chromosome addition lines with individual maize (Zea mays L.) chromosomes present in plants with a complete oat (Avena sativa L.) chromosome complement provides a unique opportunity to analyze the organization of centromeric regions of each maize chromosome. A DNA sequence, MCS1a, described previously as a maize centromere-associated sequence, was used as a probe to isolate cosmid clones from a genomic library made of DNA purified from a maize chromosome 9 addition line. Analysis of six cosmid clones containing centromeric DNA segments revealed a complex organization. The MCS1a sequence was found to comprise a portion of the long terminal repeats of a retrotransposon-like repeated element, termed CentA. Two of the six cosmid clones contained regions composed of a newly identified family of tandem repeats, termed CentC. Copies of CentA and tandem arrays of CentC are interspersed with other repetitive elements, including the previously identified maize retroelements Huck and Prem2. Fluorescence in situ hybridization revealed that CentC and CentA elements are limited to the centromeric region of each maize chromosome. The retroelements Huck and Prem2 are dispersed along all maize chromosomes, although Huck elements are present in an increased concentration around centromeric regions. Significant variation in the size of the blocks of CentC and in the copy number of CentA elements, as well as restriction fragment length variations were detected within the centromeric region of each maize chromosome studied. The different proportions and arrangements of these elements and likely others provide each centromeric region with a unique overall structure.

Stochastic processes strongly influence HIV-1 evolution during suboptimal protease-inhibitor therapy

It has long been assumed that HIV-1 evolution is best described by deterministic evolutionary models because of the large population size. Recently, however, it was suggested that the effective population size (Ne) may be rather small, thereby allowing chance to influence evolution, a situation best described by a stochastic evolutionary model. To gain experimental evidence supporting one of the evolutionary models, we investigated whether the development of resistance to the protease inhibitor ritonavir affected the evolution of the env gene. Sequential serum samples from five patients treated with ritonavir were used for analysis of the protease gene and the V3 domain of the env gene. Multiple reverse transcription–PCR products were cloned, sequenced, and used to construct phylogenetic trees and to calculate the genetic variation and Ne. Genotypic resistance to ritonavir developed in all five patients, but each patient displayed a unique combination of mutations, indicating a stochastic element in the development of ritonavir resistance. Furthermore, development of resistance induced clear bottleneck effects in the env gene. The mean intrasample genetic variation, which ranged from 1.2% to 5.7% before treatment, decreased significantly (P < 0.025) during treatment. In agreement with these findings, Ne was estimated to be very small (500–15,000) compared with the total HIV-1 RNA copy number. This study combines three independent observations, strong population bottlenecking, small Ne, and selection of different combinations of protease-resistance mutations, all of which indicate that HIV-1 evolution is best described by a stochastic evolutionary model.

A knob-associated tandem repeat in maize capable of forming fold-back DNA segments: Are chromosome knobs megatransposons?

A class of tandemly repeated DNA sequences (TR-1) of 350-bp unit length was isolated from the knob DNA of chromosome 9 of Zea mays L. Comparative fluorescence in situ hybridization revealed that TR-1 elements are also present in cytologically detectable knobs on other maize chromosomes in different proportions relative to the previously described 180-bp repeats. At least one knob on chromosome 4 is composed predominantly of the TR-1 repeat. In addition, several small clusters of the TR-1 and 180-bp repeats have been found in different chromosomes, some not located in obvious knob heterochromatin. Variation in restriction fragment fingerprints and copy number of the TR-1 elements was found among maize lines and among maize chromosomes. TR-1 tandem arrays up to 70 kilobases in length can be interspersed with stretches of 180-bp tandem repeat arrays. DNA sequence analysis and restriction mapping of one particular stretch of tandemly arranged TR-1 units indicate that these elements may be organized in the form of fold-back DNA segments. The TR-1 repeat shares two short segments of homology with the 180-bp repeat. The longest of these segments (31 bp; 64% identity) corresponds to the conserved region among 180-bp repeats. The polymorphism and complex structure of knob DNA suggest that, similar to the fold-back DNA-containing giant transposons in Drosophila, maize knob DNA may have some properties of transposable elements.

Copy-up mutants of the plasmid RK2 replication initiation protein are defective in coupling RK2 replication origins.

The broad host range plasmid RK2 replicates and regulates its copy number in a wide range of Gram-negative bacteria. The plasmid-encoded trans-acting replication protein TrfA and the origin of replication oriV are sufficient for controlled replication of the plasmid in all Gram-negative bacteria tested. The TrfA protein binds specifically to direct repeat sequences (iterons) at the origin of replication. A replication control model, designated handcuffing or coupling, has been proposed whereby the formation of coupled TrfA-oriV complexes between plasmid molecules results in hindrance of origin activity and, consequently, a shut-down of plasmid replication under conditions of higher than normal copy number. Therefore, according to this model, the coupling activity of an initiation protein is essential for copy number control and a copy-up initiation protein mutant should have reduced ability to form coupled complexes. To test this model for plasmid RK2, two previously characterized copy-up TrfA mutations, trfA-254D and trfA-267L, were combined and the resulting copy-up double mutant TFrfA protein TrfA-254D/267L was characterized. Despite initiating runaway (uncontrolled) replication in vivo, the copy-up double-mutant TrfA protein exhibited replication kinetics similar to the wild-type protein in vitro. Purified TrfA-254D, TrfA-267L, and TrfA-254D/267L proteins were then examined for binding to the iterons and for coupling activity using an in vitro ligase-catalyzed multimerization assay. It was found that both single and double TrfA mutant proteins exhibited substantially reduced (single mutants) or barely detectable (double mutant) levels of coupling activity while not being diminished in their capacity to bind to the origin of replication. These observations provide direct evidence in support of the coupling model of replication control.

In vivo function of mutated spliced leader RNAs in Caenorhabditis elegans

The role of spliced leader RNA (SL RNA) in trans-splicing in Caenorhabditis elegans has been studied through a combination of in vitro mutagenesis and in vivo complementation of rrs-1 mutant nematodes, which lack endogenous SL1 RNA. Three classes of mutant SL1 RNAs have been found—those that rescue the lethal phenotype at low concentration of transforming DNA, those that rescue at high but not low concentration, and those that do not rescue at all. These studies showed that some mutations in the otherwise highly conserved 22-nt spliced leader are tolerated for splicing and post-splicing events. A longer spliced leader also can be tolerated but only when present in high copy number. Changes in the first 16 nucleotides result in the appearance of no SL RNA, consistent with the in vitro studies by others showing that the SL1 RNA promoter partly resides within the spliced leader sequence.

Segregation of viral plasmids depends on tethering to chromosomes and is regulated by phosphorylation

Eukaryotic viruses can maintain latency in dividing cells as extrachromosomal nuclear plasmids. Segregation and nuclear retention of DNA is, therefore, a key issue in retaining copy number. The E2 enhancer protein of the papillomaviruses is required for viral DNA replication and transcription. Viral mutants that prevent phosphorylation of the bovine papillomavirus type 1 (BPV) E2 protein are transformation-defective, despite normal viral gene expression and replication function. Cell colonies harboring such mutants show sectoring of viral DNA and are unable to maintain the episome. We find that transforming viral DNA attaches to mitotic chromosomes, in contrast to the mutant genome encoding the E2 phosphorylation mutant. Second-site suppressor mutations were uncovered in both E1 and E2 genes that allow for transformation, maintenance, and chromosomal attachment. E2 protein was also found to colocalize to mitotic chromosomes, whereas the mutant did not, suggesting a direct role for E2 in viral attachment to chromosomes. Such viral hitch-hiking onto cellular chromosomes is likely to provide a general mechanism for maintaining nuclear plasmids.

Genetic analysis using genomic representations

Analysis of the genetic changes in human tumors is often problematical because of the presence of normal stroma and the limited availability of pure tumor DNA. However, large amounts of highly reproducible “representations” of tumor and normal genomes can be made by PCR from nanogram amounts of restriction endonuclease cleaved DNA that has been ligated to oligonucleotide adaptors. We show here that representations are useful for many types of genetic analyses, including measuring relative gene copy number, loss of heterozygosity, and comparative genomic hybridization. Representations may be prepared even from sorted nuclei from fixed and archived tumor biopsies.

Phenotypic conversion of drug-resistant bacteria to drug sensitivity

Plasmids that contain synthetic genes coding for small oligoribonucleotides called external guide sequences (EGSs) have been introduced into strains of Escherichia coli harboring antibiotic resistance genes. The EGSs direct RNase P to cleave the mRNAs transcribed from these genes thereby converting the phenotype of drug-resistant cells to drug sensitivity. Increasing the EGS-to-target mRNA ratio by changing gene copy number or the number of EGSs complementary to different target sites enhances the efficiency of the conversion process. We demonstrate a general method for the efficient phenotypic conversion of drug-resistant bacterial cultures.